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Abstract 

The current early stage in the investigation of the stability of the Kerr metric is 
characterized by the study of appropriate model problems. Particularly interesting is 
the problem of the stability of the solutions of the Klein-Gordon equation, describing 
the propagation of a scalar field in the background of a rotating (Kerr-) black hole. 
Results suggest that the stability of the field depends crucially on its mass /i. Among 
others, the paper provides an improved bound for fi above which the solutions of the 
reduced, by separation in the azimuth angle in Boyer-Lindquist coordinates, Klein- 
Gordon equation are stable. Finally, it gives new formulations of the reduced equation, 
in particular, in form of a time-dependent wave equation that is governed by a family 
of unitarily equivalent positive self-adjoint operators. The latter formulation might 
turn out useful for further investigation. On the other hand, it is proved that from the 
abstract properties of this family alone it cannot be concluded that the corresponding 
solutions are stable. 

1 Introduction 



Kerr space-time is the only possible vacuum exterior solution of Einstein's field equations 
describing a stationary, rotating, uncharged black hole with non-degenerate event horizon 
[31] and is expected to be the unique, stationary, asymptotically flat, vacuum space-time 
containing a non-degenerate Killing horizon [2]. Also, it is expected to be the asymptotic 
Umit of the evolution of asymptotically flat vacuum data in general relativity. 
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An important step towards establishing the vahdity of these expectations is the proof of the 
stability of Kerr space-time. In comparison to Schwarzschild space-time, where linearized 
stability has been proved, this problem is complicated by a lower dimensional symmetry 
group and the absence of a Killing field that is everywhere time-like outside the horizon. 
For instance, the latter is reflected in the fact that energy densities corresponding to the 
Klein-Gordon field in a Kerr gravitational field have no definite sign. This absence com- 
plicates the application of methods from operator theory and of so called "energy methods" 
that are both employed in estimating the decay of solutions of hyperbolic partial differen- 
tial equations.' 

On the other hand, two facts are worth noting. For this, note that in the following any 
reference to coordinates implicitly assumes use of Boyer-Lindquist coordinates [8]. 



First, in addition to its Killing vector fields that generate one-parameter groups of sym- 
metries (isometrics), Kerr space-time admits a Killing tensor [33] that is unrelated to its 
symmetries. Initiated by his groundbreaking work [10] on the complete separability of 
the Hamilton-Jacobi equation in a Kerr background. Carter discovered that an operator 
that is induced by this Killing tensor commutes with the wave operator On the other 
hand. Carter's operator contains a second order time derivative [11]. An analogous op- 
erator has been found for the operator governing linearized gravitational perturbations of 
the Kerr geometry [20]. A recent study finds another such 'symmetry operator' which 
only contains a first order time derivative and commutes with a rescaled wave operator 
[7]. Differently to Carter's operator, this operator is analogous to symmetry operators 
induced by one-parameter group of isometrics of the metric, in that it induces a map- 
ping in the data space that is compatible with time evolution, and therefore describes a 
true symmetry of the solutions. It is likely that an analogous operator can be found for a 
rescaling of the linearized operator governing gravitational perturbations of the Kerr ge- 
ometry. In case of existence, it should facilitate the generalization to a Kerr background of 
the Regge-Wheeler-Zerilli-Moncrief (RWZM) decomposition of fields on a Schwarzschild 
background [30, 35, 26, 32, 28, 15] which in turn should gready simpHfy the analysis of 
the stability of Kerr space-time. 



Second, there is a Killing field that is time-like in an open neighborhood of the event 
horizon given by 

^=a. + ^a., (1.0.1) 

where dt , d^p are coordinate vector fields of Boyer-Lindquist coordinates corresponding to 
the coordinate time t and the azimuthal angular coordinate (p, M > is the mass of the 



For the first, see, for instance, [5]. For the second, see, for instance. Chapter 2 of [27]. 
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black hole and a e [0, M] its rotational parameter. Moreover, if 

^ is time-like in the ergoregion, see Lemma 2.2. On the other hand, dt itself is space- 
like in the ergoregion, null on the stationary limit surface and time-like outside. For these 
reasons, at least for a satisfying (1.0.2), it might be possible to "join" energy inequalities 
belonging to the Killing fields by ^ and dt ■ 

The discussion of the stability of the Kerr black hole is in its early stages. The first in- 
termediate goal is the proof or disproof of its stability under "small" perturbations. As 
mentioned before, the linearized stability of the Schwarzschild metric has already been 
proved. In that case, by using the RWZM decomposition of fields in a Schwarzschild 
background, the question of the stability can be completely reduced to the question of 
the stability of the solutions of the wave equation on Schwarzschild space-time. For Kerr 
space-time, a similar reduction is not known. If such reduction exists, there is no guarantee 
that the relevant equation is the scalar wave equation. It is quite possible that such equa- 
tion contains an additional (even positive) potential term that, similar to the potential term 
introduced by a mass of the field, could result in instability of the solutions. Second, an 
instability of a massive scalar field in a Kerr background could indicate instability of the 
metric against perturbations by matter which generically has mass. If this were the case, 
even a proof of the stability of Kerr space-time could turn out as a purely mathematical 
exercise with little relevance for general relativity. Currently, the main focus is the study of 
the stability of the solutions of the Klein-Gordon field on a Kerr background with the hope 
that the results lead to insight into the problem of linearized stability. Although the results 
of this paper also apply to the case that /i = 0, its main focus is the case of Klein-Gordon 
fields of mass /i > 0. 

Quite differently from the case of a Schwarzschild background, the results for these test 
cases suggest an asymmetry between the cases ^ = Q and /i 7^ 0. In the case of the wave 
equation, i.e., ji ~ results point to the stability of the solutions [34, 16, 12, 1, 23, 24], 
whereas for /i 0, there are a number of results pointing in the direction of instability of 
the solutions under certain conditions [13, 14, 36, 17, 22, 9, 19]. 

In particular, unstable modes were found by the numerical investigations by Furuhashi 
and Nambu for ^iM ^ 1 and (a/M) = 0.98, by Strafuss and Khanna for fiM ~ 1 and 
{a/M) ^ 0.9999 and by Cardoso and Yoshida for fiM ^ 1 and 0.98 ^ (a/M) < 1. The 
analytical study by Hod and Hod finds unstable modes for fiM ^ 1 with a growth rate 
which is four orders of magnitude larger than previous estimates. On the other hand, [3] 
proves that the restrictions of the solutions of the separated, in the azimuthal coordinate. 



3 



Klein-Gordon field (RKG) are stable for 



2Mr+ 



\m\a 




(1.0.3) 



Here to e Z is the 'azimuthal separation parameter' and := M + ^/AP — a?. So far, 
this has been the only mathematically rigorous result on the stability of the solutions of the 
RKG for /^t > 0. This result contradicts the result of Zouros and Eardley, but is consistent 
with the other results above. In addition, there is the numerical result by Konoplya and 
Zhidenko, [25] which confirms the result of Beyer, hut also finds no unstable modes of the 
RKG for fiM < 1 and fiM - 1. 

Among others, this paper improves the estimate (1.0.3). It is proved that the solutions 
of the RKG are stable for fi satisfying 



Further, it gives new formulations for RKG, in particular, in form of a time-dependent 
wave equation that is governed by a family of unitarily equivalent positive self-adjoint op- 
erators. The latter might turn out useful in future investigations. On the other hand, it is 
proved that from the abstract properties of this family alone it cannot be concluded that 
the corresponding solutions are stable. 

The remainder of the paper is organized as follows. Section 2 gives the geometrical setting 
of the discussion of the solutions of the RKG and a proof of the above mentioned property 
of the Killing field ^. Section 3 gives basic properties of operators read off from the equa- 
tion, including some new results. These properties provide the basis for a formulation of 
the initial-value problem for the equation in Section 4 which is less dependent on methods 
from semigroups of operators than that of [3]. Section 4 also contains the improved result 
on the stability of the solutions of RKG, a formulation of the RKG in terms of a time- 
dependent wave equation and the above mentioned counterexample. Finally, the paper 
concludes with a discussion of the results and 2 appendices that contain proof of results 
that were omitted in the main text to improve the readability of the paper 

2 The Geometrical Setting 

In Boyer-Lindquist coordinates', {t, r,d,ip) : i7 — > R'*, the Kerr metric g is given by 

9 = 9tt dt^dt + gtip {dt ® dip + dip® dt) + grr dr (g) dr + gee dO ® dO + g^^ dip (g) d(p , 

' If not otherwise indicated, tlie symbols t, r, 9, <p denote coordinate projections wliose domains will be obvious 
from the context. In addtion, we assume the composition of maps, which includes addition, multiplication 
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where 



9u ■■= 1 



2Mr 
AS . 



2 Mar sin^ 



1 9r r ^ 



:= -s 



■ sm 



M is the mass of the black hole, a £ [0, M] is the rotational parameter and 



A := - 2Mr + , E := + cos^i 

A A A 



r+ := M + ^JaP - a2 , r_ := i\/ - ^/Ip 



X (r+,oo) X (0,7r) x (— 7r,7r) . 



In these coordinates, the reduced Klein-Gordon equation corresponding to m e Z, gov- 
erning solutions i/j : — > C of the form 



where u : fi. 



Qs := (r+,oo) x (0,7r) 
for all t E M., If E (— tt, tt), (r, 9) E ilg, is given by 



where 



^ 4mAfar 
■~ AE 



5^ 



(2.0.4) 



4mAf ar 



4TOMar 



(r2 + a2)2 - a2 A sin^ 9 (r^ + a2)E + 2Ma2r sin^6' 



E \ or or 



9 ? 

ma 

A 



1 g . nd_ 

sin 9 89 89 



sm 



for every / G C^(r2s, C) and yu > is the mass of the field. In particular, note that h 
defines a real-valued bounded function on fig which positive for m > and negative for 
m < 0. For this reason, it induces a bounded self-adjoint (maximal multiplication) opera- 
tor B on the weighted L^-space X, see below, which is positive for m > and negative 



and so forth, always to be maximally defined. For instance, the sum of two complex-valued maps is defined 
on the intersection of their domains. Finally, we use Planck units where the reduced Planck constant h, the 
speed of light in vacuum c, and the gravitational constant 7, all have the numerical value 1. 
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for m < 0. Further, D^g is singular since the continuous extensions of the coeffcients of 
its highest (second) order radial derivative vanish on the horizon {r+} x [0, tt]. 

In particular, the following proves that the Killing field 



2Mr4 



is time-like in an open neighborhood of the event horizon and time-like in the ergoregion 
if 

Af 3 ' 

Proofs are given in Appendix 1. 

Lemma 2.1. Let M > 0, a > 0. For every s £ R, the function 

g{dt + sd^,dt + sd^) 
has a continuous extension to fig. This extension is positive on 90s if and only if 



2Mr4 



Further, 



is time-like precisely on 



2Mr+ 



2Mr+ - c? €iY?e - a A^/^ sin 6i ( 1 





2M Y 








r — T- J 



e2 

Proof. See Appendix 1. 

Lemma 2.2. Let M > 0, a > and fiei, defined by 

VL^i := (a^ sin^e* - A)-i((0, oo)) , 



((0,oo)) 



□ 



denote the ergoregion. If 



then 



M - 3 



(2.0.5) 



Proof. See Appendix 1. 



□ 
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3 Basic Properties of Operators in the Equation 

In a first step, we represent (2.0.4) as a differential equation for an unknown function u 
with values in a Hilbert space. For this reason, we represent formal operators present 
in (2.0.4) as operators with well-defined domains in an appropriate Hilbert space and, 
subsequently, study basic properties of the resulting operators. Theorems 3.5, 3.6 provide 
new results. 

Definition 3.1. In the following, X denotes the weighted i^-space X defined by 

X := Ll{n,,^sme) . (3.0.6) 

Further, B is the bounded linear self-adjoint operator on X given by 

B.f := bf (3.0.7) 

for every f E X. Note that B is positive for m > and negative for m < 0. 

Remark 3.2. We note that, as consequence of the fact that B e L{X, X) is self-adjoint, 
the operator 

exp((iV2)B) , 

where exp denotes the exponential function on L{X,X), see, e.g.. Section 3.3 in [5], is 
unitary for every t e M and coincides with the maximal multiplication operator by the 
function exp{{it/2)b). 

Definition 3.3. (Definition of ^o) 

(i) We define D{Aq) to consist of all / G C^(f2s, C) H X satisfying the conditions a), 
b) and c): 

a) D^J e X, 

b) there is i? > such that /(r, 0)^0 for all r > i? and 6 e le '-^ (0, tt), 
c) 



lim ^{r,e)=0 



for all 61 e Ig. 
(ii) For every / e D{Ao), we define 

Aof Dl,f . 

Lemma 3.4. Aq is a densely-defined, linear, symmetric and essentially self-adjoint oper- 
ator in X. In addition, the closure Aq of Aq is semibounded with lower bound 
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Proof. See Lemma 2 and Theorem 4 in [7]. □ 
Theorem 3.5. The span, D, of all products 

/ ® {pr o cos) , 

where / G Cq ((r+, oo), C) and P™ : (— 1, 1) — > R is the generalized Legendre polyno- 
mial corresponding to to € Z and I e { |to| , |?t!,| + 1, . . . }, is a core for Aq. 

Proof. For this, we use the notation of [7]. According to the proof of Theorem 4 of [7], the 
underlying sets of X and X i^(f2s, (r^/A) sin 6*)) are equal; and the norms induced 
on the common set are equivalent, the maximal multiplication operator T^4/(^^-^-j by the 

function r^/(AS]) is a bijective bounded linear operator on X that has a bounded linear 
inverse; the operator H, related to Aq by 

Aq ~ T^r-^i{i^)^ I (3.0.8) 

is a densely-defined, linear, symmetric, semibounded and essentially self-adjoint operator 
in X, and D is contained in the (coinciding) domains of Aq and H. Further, it has been 
shown that (H — \)D is dense in X for A < /3, where /3 := — TO^a^/r^ is a lower bound 
for H. From this follows that D is a core for the closure H of H. For the proof, let 
/ G D{H). Since [H — \)D is dense in X, there is a sequence /i, /2, . . . in I? such that 

lim {H-X)f, = {H-\)f . 

Since H — Xis bijective with a bounded inverse, the latter implies that /i, /2, . . . is con- 
vergent to / and also that 

lim Hf„ = Hf . 

Hence, we conclude that H coincides with the closure of iJ I n. Since T,4 / / a v> , T~^, , ^— , G 

I r j ^iAz^j (AS) 

X), from the latter also follows that Aq coincides with the closure of |d- CH 
Theorem 3.6. The operator Aq coincides with the Friedrichs extension of the restriction 

ofAotoCff (as,c). 

Proof. As a consequence of Theorem 3 in [3], it follows that D is contained in the domain 
of the Friedrichs extension Ap of the restriction of Aq to C|^(f2s, C) and that ^4^?/ ^ AqJ 
for every / G I?. In this connection, note that the addition of a multiple of the identity 
operator 'does not affect' the Friedrichs extension of an operator ' Since D is a core for 
Aq, from this follows that Ap ^ Aq and hence, since yip is in particular symmetric and 
Aq is self-adjoint, that Ap = Aq. □ 



' I.e., if A is a densely-defined, linear, symmetric and semibounded operator in some Hilbert space X and 
7 G M, then the Friedrichs extension of A + 7, (A + 7)^?, and the sum of the Friedrichs extension of A, Ap, 
and 7 coincide, (A + 7)f = + 7- 
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Lemma 3.7. 

A:=Ao + {l/A) 
is a densely-defined, linear and positive self-adjoint operator in X. 

Proof. That ^ is a densely-defined, linear and self-adjoint operator in X is a consequence 
of Theorem 3.4 and the Rellich-Kato theorem. For the latter, see e.g. Theorem X.12 in 
[29], Vol. II. The positivity of ^ is a simple consequence of the fact that 



A 



sm 



^ u2 2 

- = m 
4 



A - sin^ 



AS sin^ 9 



(As; 



( AE )2 sin^ e 



[(A - sin^ 9) AS + Ul^a^r^ sin^ 



( AS )2 sin-' 9 



{ (A - sin^ 9) [A(S + 2Mr) + AM^r^] + AM^a^r^ 5111^9 } 



[ (A - sin^ 9) (S + 2Mr) + iM^r^] 



AS sin^6l 

™' [ (S - 2A/r) (S + 2Mr) + 4M\^] = 



AS^ sin^ 61 



> 



AS sin^6' 



□ 



4 Formulation of an Initial Value Problem 

In the following, we give an initial value formulation for equations of the type of (2.0.4) 
whose possibility is indicated by Theorem 4.11 in [4], see also Theorem 5.4.11 in [5]. 
Here, we give the details of such formulation, including abstract energy estimates that 
provide an independent basis for the estimate (1.0.3) and also for its improvement (4.0. 1 3) 
below. Specialization of the abstract formulation to X given by (3.0.6), A := Aq ~ C, B 
given by (3.0.7) and C :— — (a + e) for some e > 0, provides an initial-value formulation 
for (2.0.4) on every open interval / of R along with quantities that are conserved under 
time evolution. Note that in this case A + C = Aq. For convenience, the proofs of the 
following statements are given in the Appendix 2. 

Assumption 4.1. In the following, let {X, ( | ) ) be a non-trivial complex Hilbert space and 
A be a densely-defined, linear and strictly positive self-adjoint operator in X. 

Definition 4.2. We denote by W\ the complex Hilbert space' given by D{A^^^) equipped 
with the scalar product ( | ) j^, defined by 

' W\ may be regarded as a generalized Sobolev space. 
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for every ^, 7y <E D{A^^'^), and induced norm j| || i. 
Remark 4.3. Note that, as a consequence of 

for every ^ € D{A^/^), the imbedding W\ ^ X is continuous. 

Assumption 4.4. Let B : £)(A^/^) X be a symmetric linear operator in X for which 
there are a G [0, 1) and b G [0, oo) such that 

for every ^ € D{A'-/'^). Note that this imphes that B € L(WX,X). Further, let C G 
L{W\, X) be a symmetric linear operator in X and / be a non-empty open interval of K. 

Deflnition 4.5. We define a solution space Sj to consist of all differentiable u : I ^ W\ 
with Ran(u) C D{A), such that u' : I ^ X is, differentiable and 

{u')'{t) + iBu'{t) + {A + C)u{t) = (4.0.9) 

for every t € 1} 

Note that (4.0.9) contains two types of derivatives. Every first derivative of u is to be under- 
stood in the sense of derivatives of -valued functions, whereas every further derivative 
is to be understood in the sense of derivatives of X-valued functions. Unless otherwise 
indicated, this convention is also adopted in the subsequent part of this section. On the 
other hand, since the imbedding W\ '-^ X m continuous, differentiability in the sense 
of lyi^ -valued functions also implies differentiability in the sense of X-valued functions, 
including agreement of the corresponding derivatives. In particular, every u E Sj also 
satisfies the equation 

u"{t) + iBu'{t) + {A + C)u{t) = (4.0.10) 

for every t £ I, where here all derivatives are to be understood in the sense of derivatives 
of X-valued functions. Further, note that the assumptions on C, in general, do not imply 
that A + C is self-adjoint. 

Remark 4.6. According to Theorem 4.11 in [4], see also Theorem 5.4.11 in [5], for every 
to E I, £, E D{A) and 77 e W\, there is a uniquely determined corresponding u E Sj such 
that u{to) = ^ and u'{to) = 1]. The proof uses methods from the theory of semigroups of 
operators. Independently, the uniqueness of such u follows more elementary from energy 
estimates in part (iii) of the subsequent Lemma 4.7. 

' Note that the differentiabiUty of u unpUes that Ranu' C W^. 
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Parts (i) and (ii) of the subsequent Lemma 4.7 give a "conserved current" and a "conserved 
energy", respectively, that are associated with solutions of (4.0.9). Part (iii) gives associ- 
ated energy estimates, that, in particular, imply the uniqueness of the initial value problem 
for (4.0.9) stated in (iv). 

Lemma 4.7. Let u e Sj and to e /. Then the following holds. 

(i) If w g Sj, then j^ y : / — > C, defined by 

JuAt) {uit)\v'it)) - {u'{t)\v{t))+t{uit)\Bv{t)) 
for every t G /, is constant. 

(ii) The function R, defined by 

E^{t) := \\u'{t)f + {u{t)\{A + C)u{t)) 
for every t e /, is constant. 

(iii) In addition, let A + C be semibounded with lower bound 7 e R. Then 



Mt2)\\ < { 



[\\uih)\\ + _t,)]el7l^^^(t2-tl) if ^ < 

\\uih)\\+E^/^{t2-h) if7 = 

(2S„/7)i/2 (^i - e-7^/"(t2-ti)^ + ||M(ti)||e-^'^'(*=-*i) if 7 > 



for ti,t2 E I such that ti < t2- 
(iv) In addition, let A + C be semibounded. If v E Si is such that 

u(to) = w(to) , u'lta) = v'(to) , 

then V = u. 

Proof. See Appendix 2. □ 

The following example proves that it is possible that the energy assumes strictly negative 
values, but that the solutions of (4.0.9) are stable, i.e., that there are no exponentially 
growing solutions. This is different from the case of vanishing B, where there are unstable 
solutions of (4.0.9) if and only if the energy assumes strictly negative values. 

Example 4.8. The example uses for the Hilbert space X the space equipped with 
the Euclidean scalar product, A := A + C and B are the linear operators on whose 
representations with respect to the canonical basis are given by the matrices 

^ and I) , (4.0.11) 



-1 — VI 3 
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respectively. An analysis shows that A and B are bounded linear and self-adjoint operators 
in X, A is semibounded, B is positive and A+ {1/A)B'^ is strictly positive. Further, A and 
B do not commute. Finally, the operator polynomial (C — > L{X, X), A — XB — A^) 
has 4 distinct real eigenvalues. Therefore, in this case, there are no exponentially growing 
solutions of the corresponding equation (4.0.9). Fig 1 gives the graph of p := (M — > 
L{X, X),Xi^ det{A - AS - A^)) = A"* + 6A^ + SA^ - 1 which suggests that there are 
precisely 4 distinct real roots. Indeed, we notice that 

pi-5) > , p(-4) < , p{-l) > , p{0) < , p{l) > 

and hence that p has real roots in the intervals (—5,-4), (—4,-1), (—1,0) and (0,1). In 
addition, the value of the conserved energy £'„ corrresponding to the solution u of (4.0.9) 
with initial datau(O) = *(0, 1) andw'(O) *(0,0) is < 0. 

There are other possible definitions for the energy that is associated with solutions of 
(4.0.9). In cases of vanishing B, such are usually not of further use. In the case of a 
nonvanishing B, they can be useful as is the case for the RKG. In this case, the positivity 
of Es^u for sufficiently large masses of the field and 

(4.0.12) 



2Mr+ 



provides a basis for (1.0.3) and its improvement (4.0.13) below. 
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Corollary 4.9. Let s e M and u e 5/. Then, the function E^^u : / M, defined by 

Es,u{t) — \W'{t) + + {u{t)\{A + C + s{B- s))u{t)) 

for every i G /, is constant. If A + C + s(i? — s) is additionally semibounded with lower 
bound 7 G M, then 

r [ \\u{t,)\\ + |i?,,„|i/2(i2 - h) ]el^l''' if 7 < , 

\\u{t2)\\ < I \\u{h)\\ + eI!^ {t2 ^ h) if7-0 , 

[(2i?,,j7)V2 (i _ e-'''''^*^-*^)) + ||w(ti)||e-^'''(*^-*i) if 7 > , 

for ti,t2 E I such that ii < 

Proof. See Appendix 2. □ 

Theorem 4.10. If there is s e M such that A + C + s{B — s) is positive, then there are 
no exponentially growing solutions of (4.0.9). 

Proof. The statement is a direct consequence of Corollary 4.9 (or Theorem 4.17 (ii) in [4], 
see also Theorem 5.4.17 (ii) in [5]). 

□ 

Assumption 4.11. In the following, we assume that X is given by (3.0.6), A := Aq — C, 
B is given by (3.0.7) and C := — (a + e) for some e > 0. 

Theorem 4.10 leads to an improvement of the estimate (1.0.3). 

Theorem 4.12. If 

\m\a I 2M 
then there are no exponentially growing solutions of (4.0.9). 

Proof. Let s G M. In the following, we estimate Aq + sB — s^. For this, let / G D{Aq). 
Then 

(^o + .B-.^)/ = i(4Af -J-^sin^l + ^ + F.)/, 
\ or or sm OO oO sm 9 I 



where 



nri^a^ AmMar 
Vs := — + Ai'S + s ^ s2 E 

= + {fi^ - s2)S - 2s^Mr 
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First, we note that 



sin^ ( 



In the following, we assume that s = ma/ (2A/r+). Then 

^ _ {2sMr-maf _ /may f ma\^ 2'/lP~~^ ^ ^ 
'^^ A \ r+ / \ / r ~ T- ~ 

Further, we define 

Vs2 (m^ - s^)T. - 2s^Mr = - s^y"^ - 2s^Mr + a^{^? - s^) cof? 9 
If Ai > |s| • [1 + (2A//r+)]i/2,then 
2M 

Vs2 >s^ - 2s^Mr + a^U^ - s^) cos^ B > a^U'^ - s^) cos^ 6* > . 

As a consequence, 

1 / m2 \ 
= — ^ + > . 

Further, we conclude that 

(/ ® (Pr ° cos) I (Aa + sB- s^){f ® (P™ o cos))) 

> / siii0(/®(Pr°cos))* f-|-A|--^^sin0^ 
Jq^ \ or or smO oB oB 

if ® ° cos)) drdB > 

for every / e Cq ((r+, oo), C) and Z e {|to|, |to| + 1, . . . }. Since £> is a core for Aq, this 
implies that 

Ao + sB - > . 

Hence the statement follows from Theorem 4. 10. □ 



The following gives a connection of the operator + sB — s^, s € M, and the Killing 
field dt + sdtp. The corresponding proof is given in Appendix 2. This connection sheds 
light on the previous proof of the positivity of + sB — for s = ma/{2Mr+) for 
sufficiently large fi. Differently to gtt, the term g{dt + sdtp,dt + sd^p) is positive in a 
neighbourhood of the event horizon, but gradually turns negative away from the horizon. 
The latter is compensated by the mass term ji^p for sufficiently large fi. 

Lemma 4.13. Let s € M and ^ := + sdip. Then 

[Ao+msB - {msf]f 
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1 



5" 

for every / e -0(740), where 

p ■= -[gttg^v - [gt^f] = A 8111^6' . 

Proof. See Appendix 2. □ 

Subsequently, we rewrite (4.0.10) into an equivalent time-dependent wave equation that 
is governed by a family of unitarily equivalent positive self-adjoint operators. The latter 
equation might turn out useful for further investigation since only self-adjoint operators are 
involved. On the other hand, a subsequent example proves that from the abstract properties 
of this family alone it cannot be concluded that the solutions of the equation are stable. 

Lemma 4.14. Let B be additionally bounded and u E Si. Then, v : I ^ X defined by 

v{t) exp((it/2)B)w(t) 

for every t E I is twice differentiable in the sense of derivatives of X-valued functions 
and satisfies 

v"{t) + A{t)v{t) ^0 (4.0.14) 

for every t E I, where 

A{t) := exp{{it/2)B) (^A + C +^B^^ exp{-{it/2)B) (4.0.15) 
for every t G M. 

Proof. See Appendix 2. □ 

The previous can be used to prove the stability of the solutions of (4.0.9) in particular cases 
where the operators A + C and B commute. Note that in these cases, there is a further 
conserved "energy" associated to the solutions of (4.0.9). 

Theorem 4.15. If, in addition, ^ + C is self-adjoint and semibounded, B is bounded, 
A + C and B commute, i.e., 

Bo{A + C)D {A + C)oB , 

and 

A+C+^B\ 

is positive, then there are no exponentially growing solutions of (4.0.9). 
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Fig. 2: Graph of (R -> L{X, X), A ^ A - AB - A^) for A and B from Example 4.16. 



Proof. The statement is a simple consequence of Lemma 4. 14 and Lemma 4.7 (iii). □ 

Coming back to the statement of Lemma 4.14, for every t G I, the corresponding A{t) 
is a densely-defined, linear and self-adjoint operator in X, see, e.g.. Lemma 7.1, in the 
Appendix. In particular, if A + C + (1/4) is positive, A{t) is positive, too. For 
instance, according to Lemma 3.7, this is true in the special case of the Klein-Gordon 
equation (2.0.4). Hence in such case it might be expected that (4.0. 14) for u G Si implies 
that ||uj| is not exponentially growing since this is the case if A{t) = A for every t £ I, 
where A is a densely-defined, Unear, positive self-adjoint operator in X. In that case, u is 
given by 

u{t) = cos((t - to)A'/'Hto) + ^'(^to) (4.0.16) 

for all to,t e I, where cos((t - to)A^^'^) and sin((t - to)A'^/'^ / A'^^^) denote the bounded 
linear operators that are associated by the functional calculus for A^^^ to the restriction of 
cos((t — io)-idR) and the restriction of the continuous extension of sin((i — to)-idB)/idB 
to [0, oo), respectively, to the spectrum of A^^^ [5]. Note that the solutions (4.0.16) are 
in particular bounded if A is strictly positive. Unfortunately, this expectation is in general 
not true. A counterexample can be found already on the level of finite dimensional Hilbert 
spaces. 
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Example 4.16. The example uses for the Hilbert space X the space equipped with 
the Euclidean scalar product, A := A + C and B are the linear operators on whose 
representations with respect to the canonical basis are given by the matrices 

respectively. An analysis shows that A and B are bounded linear and self-adjoint operators 
in X, A is semibounded, B is positive and (l/4)i?^ is even strictly positive. Further, A 
and B do not commute. Finally, the operator polynomial (C — > L{X, X), Xi-^ A — XB — 
X^) has an eigenvalue with real part < 0. Therefore, in this case, there is an exponentially 
growing solution of the corresponding equation (4.0.10) and hence also of (4.0.14). Note 
that in this case, the corresponding family of operators (4.0. 15) consists of strictly positive 
bounded self-adjoint linear operators whose spectra are bounded from below by a common 
strictly positive real number. Fig 2 gives the graph of p := {M. ^ L{X, X), X i-^ dct(^ — 
XB — A^)) — + 4.6A'^ + 4.29A^ — 1 which suggests that there are precisely two distinct 
simple roots. Indeed, this is true. The proof proceeds by a discussion of the graph of p 
using the facts that 

p(-4) > , pi-3) < , p(0) < , pil) > , 
that the zeros of p ' are given by 

(-69 - Vl329)/40 , (-69 + \/l329)/40 , 

and that 

p((-69 + \/l329)/40) < . 

Thus, (C L{X, X), X i-> dct(yl — XB — A^)) has two distinct simple real roots and a 
pair of simple complex conjugate roots. 

5 Discussion 

The mathematical investigation of the stability of Kerr space-time has started, but is still 
in the phase of the study of relevant model equations in a Kerr background. The study of 
the solutions of the Klein-Gordon equation is expected to give important insight into the 
problem. 

In the case of the wave equation, i.e., for the case of vanishing mass i^i of the scalar field, 
results point to the stability of the solutions. On the other hand, inspection of the reduced 
Klein-Gordon equation, 2.0.4, reveals that the case of /i > originates from the case 
/i = by the addition of a positive bounded potential term 
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to the equation. If there were no first order time derivative present in the equation, from 
this alone it would be easy to prove that the stability of the solutions of the wave equation 
implies the stability of the solutions of the Klein-Gordon equation for non-vanishing mass. 



Even in the presence of such a derivative, it is hard to believe that the addition of such 
term causes instability. In particular, the energy estimates in Lemma 4.7, indicate a sta- 
bilizing influence of such a term. On the other hand, so far, there is no result that would 
allow to draw such conclusion. 

The numerical results that indicate instability in the case /i 7^ make quite special as- 
sumptions on the values of the rotational parameter of the black hole that do not make 
them look very trustworthy. They could very well be numerical artefacts. Moreover, the 
numerical investigation by Konoplya et al., [25], does not find any unstable modes and 
contradicts all these investigations. Also the analytical results in this area are not accom- 
panied by error estimates and therefore ultimately inconclusive. Still, apart from [36], all 
these results are consistent with the estimate on n in [3] and the improved estimate of this 
paper, above which the solutions of the reduced, by separation in the azimuth angle in 
Boyer-Lindquist coordinates, Klein-Gordon equation are stable. 

It seems that the proof of the stability of the solutions of the wave equation in a Kerr back- 
ground will soon be established. The question of the stability of the massive scalar field 
in a Kerr background is still an open problem, with only few rigorous results available, 
and displays surprising mathematical subtlety. In particular, in this case standard tools of 
theoretical physical investigation, including numerical investigations, seem too imprecise 
for analysis. Hence a rigorous mathematical investigation, like the one performed in this 
paper, seems to be enforced. 

6 Appendix 1 

In the following, we give the proofs of the Lemmatas 2. 1 and 2.2 from Section 2. 

Proof of Lemma 2. 1 . 

Proof. For this, let s e M. Then 

g{dt + sd^,dt + sd^) ^ gu + 2s gt^ + g^^ 



1 



2Mr 



+ 4s 



Marsm^e , AS 



A 



sni' 



2 



6 



+ 



S 



[-a^ + AsMar 
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A sin^6 



[-{a - 2sMrf - s^A{A + 4Mr) + a^s^A 811126*] 
[{a - 2sMr f + s^A{A + Ulr - sin^e) ] 



E S 

Hence g{dt + s d^^ , dt + s d^^) has a positive extension to the boundary of fis if and only if 



In this case. 



2Mr+ 



(a - 2sMrY + s^A{A + AMr - sm^O) 



A{A + AMr ~ sin^ 



[4M^{r - r+f + A(A + 4Mr - sm^9)] 



a^A 



A + AMr - sin^^ 



and hence 

g[dt + sd^,dt + sd^ 
A 



AM^r\Y. 
A 

AMh^ 



AM^rl - a? sxa^O \ AM"^ + A + AMr - sin'^l 



AM^ 



A 



AM^rlY. 



{2Mr+ - - a^{r - r+) sin^^ 

(2Mr+ - sin'^e)'^ - a^A sin^e ( 1 + 



2M 



Proof of Lemma 2.2. 

Proof. For this, let (r, 6) e n^i. Then 

A{r,e) < a^sin^e 



AM 



□ 
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and 



(2Mr+ - ain^ey - a^A(r, 9) sm^O 1 



4-4/ 

a sm t 



4Mr+ + A(r,0) 1 



4Mr+ + A(r,6') - 
2A/ 



>(A(r,0))2 

= (A(r,0))2-a2^(r,0) ( 1 + 
A(r, 6') 

A(r, e) 



2M 



2M 

r — r- 

2 



2M 

r — 

2" 



2„2 



2^2 



2 / 2M 
1 + 



r — 

2 



+ AMr+{Mr+ - a^) 

4 / 2M ' " 
1+ 



> 



2 ^ 2M 
1 + 



(r+ - r_) 



+ 4:Mr+{Mr+ - 
j + Mr+{Mr+-a^) 



Hence it follows that (r, 9) e rie2 if 



4 4 



^4 



(r+ - r_)4 
The latter is the case if and only if 



4 , M(r+-r_)4 , j\j2(r+-r_)^ 



< . 



2Mr, 



1 + A 1 



(r+-r_)-i 



Further, 



2A/r4 



> 



Mr. A'/2 
i > - 



1 + Wl 



2Af 2 - a2 



Hence if 



or, equivalently, if condition (2.0.5) is satisfied, it follows that (r, 6) e fle2- 
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7 Appendix 2 

In the following, we give the omitted proofs from Sections 3 and 4. 
Proof of Lemma 4.7. 

Proof, '(i)': For this, let f G J and ft, G M such that t + h e I. Then 

Ju,v 

h 

= h^^ [{u{t + h)\v'{t + h)) - {u{t + h)\v{t + h))+i {u{t + h)\Bv{t + h)) 
-{u{t)\v'{t)) + {u\t)\v{t))-z{u{t)\Bvm 

= h-^ [{u{t + h)- u{t)\v'{t + K)) + {u{t)\v'{t + h)- v'{t)) 

- {u'{t + h)\v{t + h)- v{t)) - {u'{t + h)- u'{t)\v{t)) 

+i {u{t + h) - u{t)\Bv{t + ft)) + i {Bu(t)\v(t + ft) - v{t))] . 

Hence it follows that j„ „ is differentiable in t with derivative 

jUt) = {<mv')\t)) {{u'y{t)\v{t))+^{u\t)\Bv{t))+^{Bu{t)\v'{t)) 
= {u{t)\{v')'{t) + iBv'{t)) - {{u')'{t) + iBu'{t)\v{t)) 
= - {u{t)\{A + C)v{t)) + {{A + C)u{t)\v{t)) = . 

From the latter, we conclude that the derivative of j„„ vanishes and hence that 2u,v is a 
constant function. 

'(ii)': For this, again, let t G / and ft e M such that t + hel. Further, let A := A + C. 
Then 

Eu{t + h) - Eu{t) 



= ft-i {u'{t + h)\u'{t + ft)) + {u{t + h)\Au{t + ft)) - {u'{t)\u'{t)) - {u{t)\Au{t)) 
= ft-i [{u'{t + ft) - u'{t)\u'{t + ft)) + {u\t)\u'{t + ft) - u'{t)) 

+ {u{t + h) - u{t)\Au{t + ft)) + {u{t)\A{u{t + ft) - u{t))) 
= ft-i [{u{t + ft) - u(t)\u{t + ft)) + {u{t)\u{t + ft) - u{t)) 

+ {A^''^{u{t + ft) - u{t))\A^/^u{t + ft)) + {u{t + ft) - u{t)\Cu{t + ft)) 

+ {Au{t)\u{t + ft) - u{t)) 

Hence it follows that Eu is differentiable in t with derivative 

{{u')\t)\u\t)) + {u'mu')'{t)) + {A^/^u'{t)\A^/^u{t)) + {u\t)\Cu{t)) 
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+ {iA + C)uit)\u'{t)) 

= - (iBu'it) + {A + C)u{t)\u' (t)) - {u'{t)\iBu'{t) + {A + C)u{t)) 

+ {u'{t)\{A + C)u{t)) + {{A + C)u{t)\u'{t)) 
= - {iBu'{t)\u'{t)) - {u'{t)\iBu'{t)) = . 

From the latter, we conclude that the derivative of vanishes and hence that £"„ is a 
constant function. 

'(iii)': Since A + C is semibounded with lower bound 7 € R, 

for every ^ G D{A). Hence it follows by (ii) that 

\\u'{t)f + 7!|w(i)|P = S„ - {{umA + C)u{t)) - -i\\u{t)f) < (7.0.18) 
for every t G K. If 7 = 0, the latter implies that 

\wm<Ei/^ 

for every i e /. Hence it follows by weak integration in X, e.g., see Theorem 3.2.5 in [5], 
that 



\\u'{t)\\dt<El'\t2-h) , 



where ti,t2 E I are such that ti < t2, and hence that 

\\u{h)\\<\\u{h)\\+El/'{h-h) . 

For the weak integration, note that the inclusion of W\ into X is continuous. If 7 > 0, it 
follows from (7.0. 18) along with the parallelogram identity for elements of X that 

||e-^'''*(e^'''-'*^ii)'(t)||2 = \\u{t) + -f'/Mt)f < 2( II^^'WII' + h'/VOir ) < 
and hence that 

for t g /. Hence it follows by weak integration in X that 



U{t2)-U{ti)\\ = 


/ u'{t)dt 











(tl:t2) 



< 



■'"''■''--.u)'{t)\\dt<{2E.J^Y/^ (( 
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for all ti,t2 G I such that ti < t2- The latter implies that 
Hence 

Mh)\\ < i2Ej^y/^ [ 

If 7 < 0, it follows from (7.0.18) that 

h'Wir <i?n-7h(t)ir < \Eu\+a\\u{t)f , 
for every t E I, where a := —7 > 0. The latter implies that 

IWrn < |i?„|i/2+ai/2||^(i)|| 

for every t E I. Hence it follows by weak integration in X that 



u{h)\\ 



U{t2)-U{ti)\\ = 


/ u'{t)dt 











\u'{t)\\ dt 



< \Eu\^/Ht2 ^ h) + a'/^ / \\u{t)\\dt, 
where ti, t2 S / are such that ti < t2, and 

\\u{t2)\\ < \\u{t,)\\ + |K|i/2(t2 - h) + / \\u{t)\\ dt . 

By help of the generalized Gronwall inequality from Lemma 3.1 in [18], from the latter 
we conclude that 

\Ht2)\\<[\\u{h)\\ + \E^\^'\t2-h)]e'^'''^'^-'^^ 
for ti E I and t2 E I such that <i < t2- 

'(iv)': For this, we define w :^ v — u. Then w is an element of Si such that w{to) = 
w'(to) = 0. This implies that 

E^{t) \\w'it)\\' + {wmA + C)w{t)) 

for every t E I is constant of value 0. Hence we conclude from (iii) that w{t) = Ox for 
alH G / and therefore that v = u. 

□ 

Proof of Corollary 4.9. 
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Proof. We define w : / — > W\ by 

v{t) := e'''u{t) 

for every t E I. Then f is differentiable with Ranv C and also w' : / — > X is 

differentiable such that 

v'{t) = e'''[u'{t) + isuit)] , (v'Yit) = e*"*[(u')'(0 + 2isu'{t) - s\{t)] 

for every t E I. Further, 

(v'Yit) + i{B - 2s)v'{t) + {A + C + sB - s^)v{t) 
= e"*[{u)'{t) + 2isu'{t) - s^u{t) + i{B - 2s){u'(t) + isu(t)) 

+ {A + C + sB^ s^)u{t)] 
= e''"*[(u')'(i) + 2isii'(i) - s^u{t) + iBu'{t) - 2isu'{t) - sBu{t) + 2s^u{t) 

+ {A + C + sB - s^)u{t)] 
= e'''[iu'yit) + iBu\t) + {A + C)u{t)] = 

for every t E I. Note that {X, A,B - 2s,C + sB - s^) satisfy Assumptions 4.1, 4.4. 
Hence it follows by Lemma 4.7 that the function E.^, : / R, defined by 



for every t e /, is constant. If, in addition, A + C + s{B — s) is, semibounded with lower 
bound 7 G R, then 



Proof of Lemma 4. 13. 

Proof. First, we notice that the only non-vanishing components of {g°'^){a,b)<^{t,r,e,ip}^ ^2 



E,{t) := \\v'{t)f + {v{t)\{A + C + SB- s^)v{t)) 
= \\u\t) + isu{t)\\^ + {u{t)\{A + C + sB- s^)u{t)) 



\Ht2)\\ < 



[||^;(<l)|| + |£;„|l/2(^2_^,)]el7^^^ 
\\v{h)\\ + e'J^ ih-h) 



if 7 < 
if 7 = 

'(t2-tl) if^>0 



for ti,t2 E I such that ti < t2. 



□ 



given by 



5" = 3*^=3^* = 



2Mar 



, 9 = , 9 



,66 



1 
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Further, we notice that 



where 
Hence 



tt _ dyy tip _ 9 tip ipip _ _9tt_ 

y : 9 : 9 

P P P 



p — - [gttgw - [gt^f] = A; 



tip 



gVV 

gtt '^t^'V ^-^^l 



1 



1 



9w 
1 

7^ 



9w 



As a consequence. 



gtt 



1 



m'^gtt + /i^P 
~9w 



for every / G 13 (Aq). Finally, it follows that , 

[ Aq + msB — (tos) ]/ = Aq/ + ms 2m '—f^f — {msf' f 



9' 



drV^\9"dr + -1= deV^\9''de 
V-\9\ V^Wl 
2 9 



g^p^p gw 



1 



1 



5" 

for every / e L'(Ao). 
Proof of Lemma 4. 14 



□ 



Proo/ First, if D e L{X, X) and / : / ^ X is differentiable in t e / and /i G R* such 
that t + he I, it follows that 

i [cxp((t + h)D)fit + h) - cMtD)f{t)] = cMtD) ^ [cxp{hD)fit + h) - /(t)] 
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exp{tD) 
exp{tD) 



cxp{hD) i[/(t + h)- fit)] + ^{cxp{hD)f{t) - fit)) 
expihD) (^^[fit + h)- fit)] - f'it)^ + cxpihD)f'it) 
+ ^iexpihD)fit) - fit)) 



and hence that g :=(/—> X, s exp(sZ?)/(s)) is differentiable in t with derivative 

expitD)[f'it)+Dfit)] . 

In particular, this implies, if / is twice differentiable in t e /, that g is twice differentiable 
in t with second derivative 

expitD)[rit) + 2Df'it) + D\fit)] . 

Applying the previous auxiliary result to D = (i/2)i? proves that v is twice differentiable. 
Further, from the definition of v, it follows that 

uit) = expi-iit/2)B)vit) , 

for every t E I. Application of the auxilary results above to D = — (i/2)_B leads to 

u'it) = expi-iit/2)B)(^v'it) -^Bvit) 

u"it) = expi-iit/2)B) (^v"it) - iBv'it) - i B^vit)^ . 

Hence it follows from (4.0. 10) that 

= u" it) + iBu' it) + Auit) 

= exp(-(ii/2)B) [v"it) - iBv'it) - - B^vit) + iBv'it) - iB- B 
\ 4 2 



expiiit/2)B)Acxpi-iit/2)B)vit) 
exp(-(ii/2)B) ( v"it) + i B^vit) + expiiit/2)B)Aexpi-iit/2)B)vit) 



= expi-iit/2)B) 
where A:= A + C. 



v"it) + expi{it/2)B) [A+^B^] cxp(-(ii/2)S)t;(i) 



□ 
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In the following, we give some abstract lemmatas that are applied in the text. For the 
convenience of the reader, corresponding proofs are added. 

Lemma 7.1. Let {X, ( | ) ) be a Hilbert space over IK e {R, C}, A a densely-defined, linear 
and self-adjoint operator in X and U G L{X, X) be unitary. Then, Au := [/ o A o 
is a densely-defined, linear and self-adjoint operator in X. Further, if L* < D{A) is a core 
for A, then U{D) is a core for U o Ao U^^. Also, if A is positive, then U o Ao U^^ is 
positive, too. 

Proof. First, we note that D{U o Ao U^^) = U{D{A)). Since D{A) is dense in X, for 
^ e X, there is a sequence of ^i, ^2, • ■ • of elements of D{A) such that 

lim ^U-^£,. 
1^—^00 

Hence also 

lim Ui, = i. 

As a consequence, U o Ao U^^ is densely-defined. Also, as composition of linear maps, 
U o Ao is linear. In addition, for ^, 77 e D{A), it follows that 

{U^\Uo A o U~^Uii) = {£,\Aii) = {Ai\ii) ^{UoAo U-^U£,\Ufi) 

and hence that U o Ao U^^ is symmetric. Further, if ^ G D{{U o Ao U^^)*), then 

{{UoAo U-^yClU-n) = {^\{UoAo U-^)Uj]) = {U-^ClAi]) 

for every rj € D{A). Hence ^ e U{D{A)), and 

{U-'^\Arj) ^ {AU~'(\r^) = {UAU-'^m 

for every rj € D{A). Since U{D{A)) is dense in X, this impHes that (U o Ao J7"^)*^ = 
U AU"-^^. As a consequence. 

Hence it follows that U o Ao [/"^ is self-adjoint. Further, let D < D{A) be a core for A. 
As a consequence, for every ^ G D{A) there is a sequence ^1, ^2, ■ • ■ in -D such that 

lim = lim A^^ = A^ . 

Hence US^i, U^2, ■ ■ ■ is a sequence in U{D) such that 

lim UC^ ^ Ui , lim VAV^^Ut,^ = UAU-^U^ . 

Therefore, U{D) is a core for UAU^^. Finally, if A is positive, it follows for ^ e D{A) 
that 

{ui\{u o A o u-')uo = (umM) = mo > 

and hence also the positivity of U AU'-^. □ 
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